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Abstract 



We re-examine the interesting possibility of utilizing the Glashow resonance (GR) 

r~| ' channel 17^ + e~ — >• W~ anything to discriminate between the UHE cosmic neutrinos 

originating from p'j and pp collisions in an optically thin source of cosmic rays. We propose 

■ a general parametrization of the initial neutrino flavor composition by allowing the ratios 

. $^1 /^^+ /*^7r+ ^° slightly deviate from their conventional values. A relationship 

between the typical source parameter k = + $^1 ) / + + + ) and the 

working observable of the GR Rq = + ) at a neutrino telescope is derived, 

lO . and the numerical dependence of Rq on k is illustrated by taking account of the latest 

experimental data on three neutrino mixing angles. It is shown that a measurement of Rq 

is in principle possible to identify the pure p^ interaction (k = 1), the pure pp interaction 

{k = 0) or a mixture of both of them (0 < k < 1) at a given source of UHE cosmic 

^ I neutrinos. The event rate of the GR signal against the background is also estimated. 
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1 Introduction 



The full construction of the IceCube detector [1], a km^-scale neutrino telescope at the South 
Pole, has recently been completed. It offers a great opportunity to discover ultrahigh-energy 
(UHE) cosmic neutrinos, whose existence may hopefully allow us to pin down the origin of 
UHE cosmic rays. The reason is simply that the UHE cosmic protons originating in a cosmic 
accelerator, such as a gamma ray burst or active galactic nuclei [2], unavoidably interact with 
ambient photons or protons. Such energetic pp or p7 interactions produce a large amount 
of charged pions, from which UHE cosmic neutrinos can copiously be produced. Since UHE 
cosmic neutrinos are not deflected by the interstellar magnetic field, they can be used to locate 
the cosmic accelerators if they are observed in a terrestrial neutrino telescope. 

The p7 and pp collisions at an optically thin source of UHE cosmic rays arc usually referred 
to as the conventional production mechanism of UHE cosmic neutrinos. Charged pions are 
mainly produced via p + 7 — A+ — )■ 7r+ + n in the p^ interaction or p + p ^ tt^ + X with 
X being other particles in the pp interaction [3]. So neutrinos arise from the decay chain 
TT"^ —>■//''' + i/^ —>■ e'^ + 1/^ + 17 ^^u^ and its charge-conjugate process. In an astrophysical source 
of either p7 or pp coUisions one has the same t'a + i^a flavor distribution $f : : = 1 : 2 : 0, 
where = + with and being the fluxes of and (for a — e, /j,, r) at the 
source. This initial flavor distribution is expected to change to $J : : = 1 : 1 : 1 at 
a neutrino telescope such as the IceCube, because UHE cosmic neutrinos may oscillate many 
times on the way to the Earth and finally reach a fiavor democracy [4] if the 3x3 neutrino 
mixing matrix V satisfies the \V^^\ = \V^.i\ condition (for i = 1,2, 3) [5]. Provided such a flavor 
democracy is really measured at the IceCube detector or at a more advanced neutrino telescope 
in the future, one will be essentially convinced that the measured UHE cosmic neutrinos come 
from the p~f or pp collisions (or a mixture of both of them) in a distant cosmic accelerator. 
Then an immediate and meaningful question is whether the neutrino telescope can discriminate 
between the p7 and pp interactions at the source. 

The answer to the above question is in principle affirmative, if the and fluxes can 
separately be determined at a neutrino telescope. Unfortunately, the present IceCube detector 
is unable to distinguish between the Cherenkov light patterns arising from the interactions of 
i/g and with ice. A possible way out is to detect the UHE cosmic flux by means of the 
Glashow resonance (GR) channel I7g + e~ — > W~ — > anything [6, 7], whose cross section can be 
about two orders of magnitude larger than the cross sections of V^N interactions around the 
resonant energy E-^ ~ 6.3 PeV [8]. As pointed out by Anchordoqui et al [9], the GR may serve 
for a useful discriminator of UHE cosmic neutrinos originating from p7 and pp collisions in an 
optically thin source of cosmic rays. The main purpose of the present paper is to re-examine 
this interesting possibility by paying particular attention to the flavor content of UHE cosmic 
neutrinos and its variation from a source to a telescope. 

Our work is different from the previous attempts in this connection (e.g., Ref. [4] and 
Refs. [9] — [13]) in several aspects. First, we propose a general paramctrization of the initial 
flavor distribution of UHE cosmic neutrinos originating from pj and pp collisions by allowing 
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7^ and $^^/$^+ 7^ 1. This treatment makes sense as the assumptions $^1 = (in 
the p7 interaction) and = (in the pp interaction) may not exactly hold in a realistic 
cosmic accelerator. Second, we establish an analytical relationship between three typical source 
parameters {S^ ^ K'-/Kl^ ^ and k ^ [$7, + + $^^^ +$7^ + $^1]) 

and the working observable of the GR [Rq = $^ /[^^ + $t neutrino telescope ^. Third, 

we examine the numerical dependence of Rq on k by taking account of the latest experimental 
data on three neutrino mixing angles. Our result shows that a measurement of Rq is in principle 
possible to identify the pure p7 interaction {k — 1), the pure pp interaction = 0) or a mixture 
of both of them (0 < k < 1) at a given astrophysical source, in particular after all the neutrino 
mixing parameters of V are well determined from a variety of terrestrial neutrino oscillation 
experiments. In addition, the event rate of the GR signal against the relevant background is 
also estimated in this paper. 



2 Modified Flavor Distribution on the GR 

We have denoted the tt^ fluxes from the pj interaction as , and those from the pp interaction 
as $^±- In the conventional picture of p-f collisions one mainly considers the A- resonance 
channel p + 7 — )> A"*" — )■ n + n^, and thus $^1 = is taken as a good approximation for a 
given astrophysical source. As for the pp interaction in a cosmic accelerator, the produced tt"'", 
7r~ and 7r° mesons are expected to be in almost equal amount due to the isospin symmetry. 
Hence = is also a good approximation. In general, however, a small amount of 
7r~ mesons should be produced from the p7 interaction (e.g., from the multi-pion production 
channel p + 7 — >■ n+7r''" + n(7r^7r~) with n being a positive integer [14] ^), and a slight difference 
between and must be present for the pp interaction. So we consider a general source 
in which both p7 and pp collisions are important. To be explicit, we define three typical source 
parameters to describe the content of 7r+ and 7r~ mesons produced from p7 and pp collisions: 
5^ ^ S^^ ^ - 1 and 

- $PP I ^PP I ^Pl I ^PT • 



" TT 



In this simple parametrization the k, = 1 and k, = cases correspond to the pure p7 and pure pp 
interactions, respectively. If the value of k, is found to lie in the < k < 1 range at a neutrino 
telescope, it will imply that both p^ and pp collisions exist at the relevant astrophysical source. 

Now we look at the flavor composition of UHE cosmic neutrinos originating from p^ and pp 
collisions in an optically thin source of cosmic rays. Taking account of k, 5^^ and 5^^ deflned 



■l-Note that and T have been used in Ref. [11] to describe the fraction of UHE cosmic neutrinos produced 
from the P7 interaction and the working observable at the neutrino telescope, respectively. 

§Note that the back reaction n + 7 — >^ p + 7r~ could also produce 7r~ mesons if the optical thickness of the 
source is non- negligible, and the flux originating from the beta decays of neutrons might even dominate in 
some astrophysical sources for very specific energy ranges [12]. For simplicity, here we follow Ref. [9] and focus 
on the cases in which the afore-mentioned effects can be neglected. 
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above, we obtain the ratio of neutrino and antineutrino fluxes as follows: 



( + ) { 3 : : - : - : : } + ( + ) { : - : - : - : : } 
1 
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: : 



(2) 



fiavor distribution $g : $^ : $^ 



(for a = e. fi.r), it is straightforward to arrive at the 
1:2:0. This simple result is completely 



independent of three source parameters. That is why one has to separately measure the and 
i7g fluxes at a neutrino telescope so as to probe and at the astrophysical source. 

Thanks to the effect of neutrino oscillations, the i/^ and I^^ fluxes observed at the telescope 
are simply given by 



(3) 



where = P{iy^ 



— )■ z/^) and Pq,^ = P(z/q, — )■ z/^) stand respectively for the oscillation 



probabilities of UHE cosmic neutrinos and antineutrinos. Since the galactic distances far exceed 
the observed solar and atmospheric neutrino oscillation lengths, and P^,^ are actually 
averaged over many oscillations and thus become energy-independent: 

Pa,-Pa,-i:i\yaf\V,f) , (4) 
i 

where V^^ and Vg^ (for a, f3 = e, ii,t and i — 1,2, 3) denote the elements of the 3x3 neutrino 
mixing matrix V. For our purpose, we are mainly interested in the determination of via 
the GR channel + e" — > W~ — > anything. So we estabhsh a hnk between three source 
parameters and a working observable at the neutrino telescope: 
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(5) 



en 



{2 + 5J{l + SJ 

where P^g, Pg^ and P^^ can directly be read off from Eq. (4). After the matrix elements 
of V are determined to a sufficiently good degree of accuracy in solar, atmospheric, reactor 
and accelerator neutrino oscillation experiments, a measurement of Rq at a neutrino telescope 
will allow one to constrain the source parameters via Eq. (5). There are two special cases, 
corresponding to the pmc p7 interaction {k — 1) and the pure pp interaction {k — 0) at the 
astrophysical source of cosmic rays: 
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(6) 
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If both 5^ and 5^^ are switched off, then Eq. (5) can be simphfied to 



This result is particularly interesting in the sense that it offers an opportunity to determine k 
in a cosmic accelerator from the measurement of Pq at a neutrino telescope. 

In the standard parametrization of V [15], Pgg, P^^ and P^^ can be expressed in terms of 
three neutrino mixing angles (^^2) ^23' ^13) Dirac-type CP- violating phase 5 as follows: 



P ~ 1 - 
2 



^ sin^ 26*12 - (2 - sin^ 26*12) sin^ 6*13 , 

Pg^ ~ - sin^ 29^2 cos^ ^23 + ^ ^^12 2^23 ^13 + ^2 sin^ ^23 ~ 2 ^^^^ 2^12^ sin^ 9 

11 1 
P^^ ~ 1 — - sin^ 2^23 ^ 2 ^^^^ ^^12 cos^ ^23 ~ 2 ^™ ^^12 sin 2^23 cos^ ^23 sin ^13 cos 5 



1 

+4 



sin^ 2^12 sin^ 2^23 (2 + cos 25) - 8 sin^ ^23! sin^ ^13 , (8) 



in which the terms proportional to sin^ ^^3 ~ 0.3% and those much smaller ones have been 
omitted. A global analysis of the latest neutrino oscillation data [16] yield sin^ 9-^2 — 0-306l!lo!oi5) 
sin^ ^13 = 0.0211^^^ and sin^ ^23 = 0.421^^ at the la level ^, while the Dirac-type CP- 
violating phase 5 remains entirely unrestricted. Because the contributions of 5 to P^g, Pg^ and 
P^^ are always suppressed by small sin^i3, the 5-induced uncertainties in the calculation of Rq 
should not be significant. 

Note that a real observable of the GR channel + e~ — )■ W~ anything at a neutrino 
telescope can be the ratio of the events to the and events of charged-current interactions 
in the vicinity of the resonance E-^^ ~ M^/(2mg) ~ 6.3 PeV [10, 19]: 

^ ~ A/"CC j^mCC = "-^0 ' (9) 

where a ~ 30.5 can be obtained in an optimal case by assuming the neutrino spectrum [10] 
and considering the muon events with contained vertices [18] in a water- or ice-based detector. 
A more accurate calculation of a is certainly crucial for the IceCube detector to detect the rate 
of the GR reaction [9] . Note also that the flux of E^^ ~ 6.3 PeV might largely get absorbed 
in passing through the Earth [10]. Hence it is only feasible for a neutrino telescope to detect 
the downward-going or horizontal flux whose energy lies in the vicinity of the GR, in which 
case the atmospheric neutrino flux of the same energy is negligibly small and should not be of 
concern as an important background [10]. 

We proceed to illustrate the dependence of Rq on k, 5^^ and with the help of current 
experimental data on three neutrino mixing angles. First of all, we assume 5^^ — 5^^ — Q and 



^Note that these results are obtained by using the old reactor antineutrino fluxes [16] . If the new reactor 

antincutrino fluxes [17] arc used, the corresponding best-fit values and la ranges of sin^ ^-i^d sin^ 6^^ will 
be shifted by about +0.006 and +0.004, respectively, but the result of sin^ ^23 is essentially unchanged [16]. 
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use Eq. (7) to describe the relationship between and k,. Fig. 1 shows the allowed region 
of Rq versus < k < 1, where the la ranges of 6-^2, 613 and together with 6 G [0,27r) 
have been scanned. The central value of R^ for a given value of k, is calculated by inputting 
the best-fit values of three neutrino mixing angles (i.e., sin^ 6*12 = 0.306, sin^ 6'i3 = 0.021 and 
sin^ 023 — 0.42 [16]) and taking 5 = 0. Although the uncertainties associated with four neutrino 
mixing parameters remain rather large, we have the following quantitative observations: (1) 
the magnitude of Rq is restricted to the range 0.18 < Rq < 0.58; (2) Rq lies in the range 
0.18 < Rq < 0.31 for the pure p7 interaction (i.e., k, — 1); and (3) Rq lies in the range 0.45 < 
Rq < 0.58 for the pure pp interaction (i.e., k — 0). As the neutrino mixing parameters can be 
more and more precisely measured in the ongoing and future neutrino oscillation experiments, 
we expect that the GR will serve as a clear discriminator of UHE cosmic neutrinos originating 
from pp and p7 collisions at an astrophysical source. 

Now let us examine possible effects of 6p^ and 6pp on the relationship between Rq and k,. 
For simplicity, we only take the best-fit values of three neutrino mixing angles and assume 
(5 = in our numerical illustration. The change of Rq with respect to three source parameters 
K, and is shown in Fig. 2, where 6pj G [0, +0.2] and 6pp G [—0.2, +0.2] have been 
assumed. Note that is positive (or vanishing) by definition, while 6^^ can be either positive 
or negative (or vanishing), corresponding to an excess of the 7r~ or tt"*" events (or = ^^-) 
in the pp interaction at an astrophysical source. As in Fig. 1, the central curve of Rq varying 
with K in Fig. 2 is obtained in the assumption of Sp^ = 5pp = 0. It is straightforward to see that 

and can significantly affect Rq for a given value of k. For the pure p7 interaction with 
K — 1, a variation of 5p^ from to 0.2 results in a change of Rq by more than 30% as compared 
with its original value. As indicated by Eq. (6), it is in principle possible to determine or 
constrain the free parameter 5p^ (or 5pp) for a given source with the pure p7 (or pp) interaction 
by measuring Rq at a neutrino telescope. 

If the uncertainties from both the neutrino mixing parameters {0^2^ ^13) ^23 ^^'^^ ^^^^'^ 
the source parameters {6p^ and 6pp) are taken into account, it will be almost impossible to 
distinguish between p'-f and pp collisions even if i?Q ~ 0.4 is extracted from a neutrino telescope 
experiment. This observation implies that it does make sense for us to consider the nontrivial 
effects of 6p^ and 6pp. What we can do at present is to carefully study the yields of tt^ fiuxes 
in the realistic models of p'j and pp collisions, so as to obtain some theoretical constraints on 
5p^ and 6pp [20]. In addition, we must determine the neutrino mixing parameters as precisely 
as possible in all the terrestrial neutrino oscillation experiments. 

3 Estimate of the Event Rate and Background 

To further illustrate, let us estimate the event rate of the GR signal and the relevant back- 
ground. We assume the total fiux of UHE cosmic neutrinos and antineutrinos originating from 
an optically thin source to saturate the Waxman-Bahcall (WB) bound [21] 

^2$^^_ = 2 X 10"^ GeV cm'^ s"^ sr"^ , (10) 
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where stands for the ratio of the pion energy to the initial proton energy, and ~ 3 for a 
source evolution oc (1 + 2;)'^ with z being the redshift. We have = ^ 0.25 for J97 collisions 
or = e^^ ~ 0.6 for pp collisions. Therefore, the WB bound actually depends on whether the 
pp or p7 collision is assumed. Since there is on average one cosmic-ray neutron produced per 
proton collision, we may parametrize ^^^j; saturating the WB bound as 



^ 6 X 10-^ [(1 - k') ef + K'e^] GeV cm' 



s ^ sr ^ 



(11) 



where k' denotes the fraction of the p7 collisions. In this parametrization k' — 1 and k' — Q 
correspond to the pure p7 and pure pp interactions, respectively. Note that we have defined k 
in Eq. (1) as the fraction of the pion fluxes from the p7 collisions. The relationship between 
K and k' can be easily established: 



K e. 



pp 



pp 



(12) 



Given the total flux of neutrinos and antineutrinos in Eq. (11) and their flavor distribution 
at the source in Eq. (2), it is then possible to calculate the neutrino and antineutrino fluxes 
of different flavors at a neutrino telescope by taking account of the effect of flavor oscillations. 
We obtain 
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(13) 



where Eq. (3) has been used and = 6 x 10 GeV cm s sr" -1 (1 GeV/E^Y is defined. 
Note that the energy dependence of and fiuxes in Eq. (13) has been suppressed. 

Following Ref. [9], we estimate the event rate of the GR signal in the IceCube experiment: 



dNjdt = 66.7% X TV^ffAQ J dE,^l{E,)aQ^{E^) , 



(14) 



in which the coefficient 66.7% is the branching ratio of hadronic decays, N^f^ ^ 6 x 10^^ 
denotes the number of target electrons for an effective volume V^g ~ 2 km'^ of the IccCubc de- 
tector, AQ fa 27r is the solid angle aperture, and aQ^{E^) = ng'^M^6{2m^E^ — M^) / {Am^E^) 
is the cross section of the GR scattering. The typical GR signal is the shower events induced 
by the hadronic decays oiW~ in the resonant energy region, while the main background comes 
from the non-resonant inelastic scattering of and with nucleons in the detector. As for the 
background events, the effective number of target nucleons is approximately twice the number 
of electrons (i.e., N^fi 1.2 x 10^^) and the sohd angle aperture is AQ' ?s An. The cross sections 
of charged-current v^N and V^N interactions are well represented by the power-law forms [8]: 

0.402 



'CC 



uN 
CC 



2.69 X 10 



cm2 



E,. 



2.53 X 10-=*^ cm^ 



1 GeV; 

TP . 0.404 

1 GeV; 



(15) 
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Integrating over the resonant acceptance energy bin (lO^'^ • • • lO^'^) GeV for the IceCube tele- 
scope, we can obtain the event rate for the background 

dNjdt = K^AQ' J^^^^^^^ dE^ [$J,(^.)^cc(^.) + ^liK)^lc{K)] ■ (16) 

As usual, the signal-to-background ratio can be defined as R^^^ = {dNg/dt)/{dN^/dt), which 
measures the significance of the signal events. 

We perform a numerical calculation of the event rate of the GR signal dN^/dt and the signal- 
to-background ratio i?s/b, and examine their dependence on the source parameters (k, , 5^^) 
and the neutrino mixing parameters (^^i2) ^23) ^13? Fig. 3 shows the expected event rate 
dN^/dt versus < k < 1, where 6p^ = ^pp = ^ is assumed and the Icr ranges of 6^2, ^13 and ^23 
together with 6 G [0, 27r) have been scanned. The corresponding signal-to-background ratio in 
this case is shown in Fig. 4. For the pure pp interaction (i.e., k, = 0), we obtain dA^"g/dt ^ 3.5 
per year and R^^^ ~ 5, indicating a great discovery potential of the IceCube telescope after 
several years of data accumulation [9]. For the pure pj interaction (i.e., k = 1), however, 
the event rate is quite low: dN^/dt ~ 0.8 per year. Hence it is quite challenging for the 
IceCube detector to discover UHE cosmic neutrinos originating from an optically thin source 
with the pure p7 collisions. This observation justifies the importance of the GR channel in 
distinguishing between pp and p7 interactions. On the other hand, the dependence of dN^/dt 
and Rg^y^ on the source parameters {k, 5p^, 5pp) are illustrated in Fig. 5 and Fig. 6, where 
the best-fit values of three neutrino mixing angles (i.e., sin = 0.306, sin^^is = 0.021 and 
sin^ ^23 = 0.42 [16]) together with 5 = have been input. It is straightforward to see the 
degeneracy between the uncertainty induced by those neutrino mixing parameters and that 
by the source parameters (5^^, 6pp). So a full determination of the latter requires more precise 
values of neutrino oscillation parameters and a neutrino telescope whose scale should be much 
larger than the IceCube detector. 



4 Summary 

We have re-examined the possibility of using the GR channel z7g -|- e~ — )■ W~ — )■ anything to 
discriminate between the UHE cosmic neutrinos originating from p^ and pp collisions in an 
optically thin source of cosmic rays. After proposing a general parametrization of the initial 
neutrino flavor distribution by taking account of non-zero and Spp at the source, we have 
established an analytical relationship between the typical source parameter k and the working 
observable of the GR Rq at a neutrino telescope. We have also illustrated the numerical 
dependence of Rq on k with the help of the latest experimental data on three neutrino mixing 
angles. We find that a measurement of Rq is in principle possible to identify the pure p7 
interaction {k = 1), the pure pp interaction (k = 0) or a mixture of both of them (0 < k < 1) 
at a given source of UHE cosmic neutrinos. In addition, the event rate of the GR signal against 
the relevant background is estimated by assuming the total flux of UHE cosmic neutrinos and 
antineutrinos originating from an optically thin source to saturate the WB bound. 



8 



A measurement of the GR and a determination of the flavor distribution of UHE cosmic 
neutrinos at an astrophysical source are certainly big challenges to the IceCube detector and 
other possible neutrino telescopes. Anyway, our present understanding of the production mech- 
anism of UHE cosmic neutrinos depends on a number of hypotheses and thus needs more and 
more observational supports. We therefore expect that neutrino telescopes can help us in this 
connection in the long run. 

We would hke to thank S. Pakvasa and W. Winter for their useful comments and discussions. 
This work was supported in part by the National Natural Science Foundation of China under 
grant No. 10875131 (Z.Z.X.) and by the Alexander von Humboldt Foundation (S.Z.). 
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Figure 1: The dependence of the working observable -Rq on the source parameter k in the 



assumption of 5. 



0. The dashed curve corresponds to the best-fit values of 



125 



9^^ and together with 5 = 0, and the uncertainties come from the la error bars of three 
neutrino mixing angles and an arbitrary change of 5 G [0, 2ti). 




Figure 2: The dependence of the working observable Rq on the source parameter k, where 
the best-fit values of 6^21 ^13 ^'^d ^23 together with 5 = have been taken. The dashed curve 
corresponds to 5p^ = 5. 



0, and the uncertainties come from the variations of 6^^ and 6^^ in 



the ranges G [0, +0.2] and 6pp G [-0.2, +0.2]. 
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Figure 3: The dependence of the event rate of the GR signal dN^/dt on the source parameter 
K in the assumption of Sp^ = 6pp = 0. The dashed curve corresponds to the best-fit values of 
$12, and together with 6 = 0, and the uncertainties come from the la error bars of three 
neutrino mixing angles and an arbitrary change of 5 G [0, 27r). 
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Figure 4: The dependence of the signal-to-background ratio R^^^ on the source parameter k in 
the assumption of 6p^ = 6pp = 0. The dashed curve corresponds to the best-fit values of 6*^2; 

and ^23 together with 6 = 0, and the uncertainties come from the la error bars of three 
neutrino mixing angles and an arbitrary change of 6 G [0, 27r). 
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Figure 5: The dependence of the event rate of the GR signal dN^/dt on the source parameter 
K, where the best-fit values of 6^2, ^13 and together with 6 = have been taken. The dashed 
curve corresponds to 6p^ = 6pp = 0, and the uncertainties come from the variations of 6^^ and 
6 in the ranges Sp^ G [0, +0.2] and 6pp G [-0.2, +0.2]. 
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Figure 6: The dependence of the signal-to-background ratio R^^^ on the source parameter n, 
where the best-fit values of 6-^^2j ^13 and ^23 together with 6 = have been taken. The dashed 
curve corresponds to 6py = 6pp = 0, and the uncertainties come from the variations of 6p^ and 
6 in the ranges Sp^ G [0, +0.2] and 6pp G [-0.2, +0.2]. 
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